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Abstract 

We present Monte Carlo simulations of liquid-crystalline material confined to a nanoscopic slit- 
pore. The simulations are carried out under isothermal conditions in a specialized isostress ensemble 
in which N fluid molecules are exposed to a compressional stress rn acting on the fluid in directions 
parallel with the substrate planes. Fluid-fluid and fluid-substrate interactions are modelled as in our 
previous work [M. Greschek et al, Soft Matter, 2010, DOI:10.1039/B924417D). We study several 
anchoring mechanisms at the solid substrate by introducing an anchoring function < g (u) < 1 
that depends on the orientation u of a fluid molecule relative to the substrate plane; g (u) "switches" 
the fluid-substrate attraction on or off. Here we focus on various heterogeneous (i.e., hybrid) 
anchoring scenarios imposing different anchoring functions at the opposite substrates. As in our 
previous study we compute the isostress heat capacity which allows us to identify states at which 
the confined fluid undergoes a structural transformation. The isotropic-nematic transformation 
turns out to be nearly independent of the specific anchoring scenario. This is because the director 
in the nematic phase assumes a direction that is a compromise between the directions enforced by 
the competing anchoring scenarios at either substrate. On the contrary, at lower compressional 
stresses molecules prealign in specific directions that depend on details of the anchoring scenario. 
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I. INTRODUCTION 



If a liquid crystal is in its nematic phase specific orientations of the director [i.e., the (unit) 
vector n with which molecules align preferentially] can be realized through specific anchor- 
ing scenarios at the solid substrate. The substrate-induced alignment can be perturbed if 
an external field is superimposed such that it is competing with the director. An example 
is the so-called Freedericksz cell in which a nematic liquid crystal is exposed to a magnetic 
field pointing in a direction other than that in which molecules are anchored at the walls of 
the cells 3]. If the external field is sufficiently strong and if the substrate separation is 
large enough there will be a region in the liquid crystal where the local director points in 
another direction than that enforced by substrate anchoring. In other words, the orienta- 
tional order in the nematic phase is perturbed by the presence of the external field. Because 
the correlation lengths of these perturbations may be substantial they can be used in sensor 
applications. This was recently demonstrated by Guzman et al. who considered a model 
for a bionanosensor based upon a confined liquid-crystalline phase The sensor consists 
of thin films of liquid- crystalline materials confined between substrate surfaces with specific 
anchoring characteristics. If particles such as proteins or viruses bind to these solid surfaces 
they give rise to long-range perturbations of the orientational order in a nematic phase that 
reports the binding of these particles at specific receptors on the solid surfaces 4|-|6(. The 
binding process perturbs the local structure of the liquid-crystal phase over sufficiently long 
distances so that optical methods can be used to detect the perturbation. 

The apparent importance of various anchoring scenarios for the nematic phase of such 
nanoconfined liquid crystals prompted us to investigate the impact of various anchoring 
scenarios on the isotropic-nematic (IN) phase transition in a model liquid crystal of prolate 
ellipsoids of revolution confined to a nanoscopic slit-pore Q. Our results indicate that 
the nature of the anchoring at the substrate has a pronounced effect on the location of 
the IN transition and on structural details associated with that transition. In view of our 
earlier results it thus seems interesting to replace in this work the original homogeneous 
[3] by heterogeneous (or "hybrid") anchoring scenarios such that the anchoring is different 
at the opposite substrates forming the slit-pore. This will cause the director field to be 
inhomogeneous such that interesting new physical phenomena are anticipated. 

That such heterogeneous or competing fluid-wall interactions can lead to exciting effects 
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and interesting new physics has already been demonstrated in the past in several contexts. 
Perhaps the best-known example is the localization-delocalization transition that occurs for 
the gas-liquid interface in fluids confined between two planar walls where one wall favors 
the liquid- and the other one favors the gas-like phase 8-lof . If the substrate is nonplanar 
such as in the bi-pyramid considered by Milchev et al. Uj competing surface fields can 
cause the spontaneous magnetization in an Ising magnet to vanish at a filling temperature 
Tf below the bulk critical temperature T c b; the spontaneous magnetization remains zero for 
all T > T f . 

In the context of liquid crystals surprisingly little attention has been paid so far to hybrid 



anchoring scenarios despite their importance for optical applications [12|, [13( . For example, 
Ziherl et al. consider a situation in which the molecules of a liquid crystal are aligned 
homeotropically at one wall and planar at the opposite such that in the nematic phase 
molecules appear to be frustrated with respect to the direction of n [2|. The authors focus 
on the impact of this frustration on fluctuation-induced forces that can be interpreted as an 

n 

analogue of the well-known Casimir forces in quantum systems [141 ] . For the same hybrid 
anchoring scenario Rodriguez-Ponce et al. employ mean-field density functional theory and 
observe a linear variation of the tilt angle between the local director field and the substrate 
normal (see also Ref . 1 151 ) as far as strong anchoring at the substrate is concerned 16j • Such a 
linear variation of the tilt-angle profile has also been reported by Steuer et al. [l7] for hybrid 
anchoring corresponding to a nanoscopic twisted nematic cell and for the same model system 
used in this work. However, to date a systematic study of the impact of hybrid anchoring 
on the IN transition is still lacking. 

Experimentally, hybrid anchoring scenarios can be realized in a number of ways. For 
example, Chung et al. use photoalignment to fabricate planar substrates where the anchoring 
at each substrate differs by an angle of tt/A 181 ] . They show that the resulting molecular 
orientation in a confined liquid crystal can be controlled by the photoirradiation time. The 
same authors show that very similar results can be obtained if the surfaces are rubbed in 
specific directions and under controlled conditions. Zappone et al. perform experiments 
using the surface forces apparatus (SFA) to measure the solvation force in nematic liquid 
crystals under hybrid anchoring conditions 19j. They set up a hybrid (homeotropic-planar) 
anchoring scenario by dipping the muscovite mica surfaces, with which the crossed cylinders 
in the SFA setup are coated, in various cationic surfactant solutions where the surfactants 
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differ only in the length of the aliphatic chains. 

In this work we will study the impact of various hybrid anchoring scenarios on properties 
of nanoconfined model liquid crystals by means of Monte Carlo simulations in a specialized 
isostress ensemble similar to our previous study |7|. Our emphasis will be on structural 
changes arising in the confined phase as the transverse compressional stress t\\ increases. 
Increasing tn eventually drives the system from an isotropic to a nematic phase. We have 
organized the remainder of our manuscript such that we will briefly describe the model 
system in Sec. [Til Our results are presented in Sec. IIHI and summarized and discussed in 
the concluding Sec. IIVI Because of the hybrid anchoring scenarios employed in this study 
biaxiality of the nematic phases becomes an issue. The biaxial order parameter is introduced 
in Appendix |A] using perturbation theory. 



II. THE MODEL SYSTEM 



A. Interaction potentials 



The model liquid crystal considered in this work is confined to a mesoscopic slit-pore 
such that we may decompose the total configurational energy into a fluid-fluid (ff ) and into 
a fluid-substrate (fs) contribution according to 



U{R, U) = U S (R, U) + U b (R, U) 



(2.1) 



where R = {ri, r-i, ■ ■ ■ , J*jv} arid U = {ui, U2, • • • , u n } are shorthand notations for the sets 
of center-of-mass coordinates and unit vectors specifying the orientations of the N liquid- 
crystalline molecules, respectively. More specifically, 



N N 



U S (R, U) = - ^ U S ( r ij> «i> «. 



(2.2) 



1=1 jyti 



where = — rj is the distance vector between the centers of mass of particles i and j 
assuming pairwise additivity of their interactions. As in our previous work Q we take ug 



as (see also Ref. 



2fl) 



u s i r.j.Uj.Uj, 
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> ij 



{1 + ^ (fy, U U Uj)} 



(2.3) 
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where r = \r\, and r = r/r. Hence, Ug is a Lennard- Jones potential where the attractive 
contribution is modified to account for the orientation dependence of the interaction between 
a pair of molecules. In eqn. (12. 3p . a denotes the "diameter" of a spherical reference molecule 
and e is the depth of the attractive well in that reference model. The anisotropy of the 
fluid-fluid interaction is accounted for by the function 

* (f y , u h Uj) = 5eiP 2 («i • Uj) + 5e 2 [P2 («i • rij) + A («j • f^)] (2.4) 

where we take s% = 0.04 and £9 — —0.08 as in our previous study As we have demon- 
strated there (see Fig. 1 of Ref. |7|) our molecules turn out to be prolate ellipsoids of revolution 
with an aspect ratio corresponding to 1.26. The specific functional form of ^> preserves the 
head-tail symmetry of our molecules, that is orientations tii and — Ui are equivalent. In 
eqn. O 

P 2 (x) = i (3x 2 - 1) (2.5) 

is the second Legendre polynomial. 

To model the fluid-substrate contribution to U we follow earlier work (3, and introduce 

2 N 

U is (R,U)= ^^^(r^i) (2.6) 

fc=l i=l 

where 



4s 1 (^i, Wi) = 2vr£:p s o- 2 
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c \ / o 



g [k] (fii 



(2.7) 

5 ± s z0 /2 y \^ ± Szo/2/ 
In eqn. (12. 7p . ^ is the ^-coordinate of the center-of-mass position of molecule i in a space- 
fixed Cartesian coordinate system. In this coordinate system we assume the lower substrate 
(k = 1) to be located at z w = — s z0 /2 whereas the upper one is located at z w = +s z o/2 
(k = 2). In other words, eqn. ( 12. 7ft assumes the substrates to be structureless such that 
at fixed molecular orientation «f s depends only on the distance of molecule i from either 
substrate along the z-axis. In the prefactor the areal density of the solid substrates p s = 2/£ 2 
where ija = v^4 is the lattice constant of a single (100) plane of a face-centred cubic lattice. 

B. Anchoring scenarios 

In eqn. ( 12. 7J , < (fij) < 1 is the so-called anchoring function. Choosing different 
functional forms the anchoring function permits to realize different, energetically favorable 



orientations of molecule % with respect to the substrate plane. Specifically, we take (ttj) 
to be given by one of the following expressions: 



9o W 
9± («) 
9\\ («) 
0x(«) 



we, 



(2.8a) 
(2.8b) 
(2.8c) 
(2.8d) 



The anchoring functions given in eqn. ( 12. 8ft have the following physical significance. As 
we rationalize elsewhere if (u) is given by go(u) molecules orient themselves in a 
homeotropic fashion at the substrates even though the substrates themselves do not discrim- 
inate any molecular orientation energetically. Employing the expression given in eqn. (I2.8bj) 
also results in homeotropic alignment which is, however, stronger than the one induced by 
go (u) because the substrates discriminate the homeotropic alignment directly. To overcome 
the homeotropic alignment already supported by the mere presence of the substrate we also 
consider the anchoring functions specified in eqn. (I2~8c1) and (l2~8dl) . In both planar 
orientation parallel with the substrate plane is energetically favored. The difference is that 
by employing eqn. (I2.8cl) any orientation of a molecule parallel with the substrate plane is 
energetically favorable whereas eqn. (12.8dj) selects those orientations where u is aligned with 



the x-axis. We shall distinguish homogeneous anchoring scenarios where ({^) = gl 2 ] 
from heterogeneous (i.e., "hybrid") ones characterized by (ui) ^ g^ (ui) (see Tabled]). 

TABLE I: Hybrid anchoring scenarios employed in this work. 





(u) 




k = 1 


k = 2 


anchoring scenario 


acronym 


g± («) 


9\\ (2) 


homeotropic-planar 


hp 


9o («) 


9\\ (2) 


nonspecific-planar 


np 


9o (u) 


5x (u) 


nonspecific-directional 


nd 
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III. RESULTS 



A. Numerical details 

As in our previous study [3] we employ Monte Carlo simulations in a specialized isostress 
ensemble in which the fluid is exposed to a constant compressional stress in directions 
parallel with the substrate planes. We employ periodic boundary conditions in the x- and 
y-directions. To generate a Markov chain of configurations in this ensemble we utilize the 
algorithm described in Sec. 4.1 of Ref. |21 



We refer to a Monte Carlo cycle as a sequence 
of N attempted displacements and rotations of sequentially selected fluid particles plus one 
attempted change of the area A of the computational cell in the x-y plane. Typical runs 
comprise 10 5 Monte Carlo cycles in regions where the fluid does not undergo an isotropic- 
nematic (IN) transformation. In the immediate vicinity of a structural transformation the 
length of a typical run was enlarged to 3 x 10 6 Monte Carlo cycles to guarantee sufficient 
statistical accuracy of the results (see Sec. IV B. of Ref. Q). 
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FIG. 3.1: Plot of Uf s (z, u) = u&J (z, u) + (z, u) as a function of position z/s z o in the lower half 
of the system (z/s z o < —0.5) for s z o = 20. The plots are symmetric with respect to the line z = 
and were obtained assuming g 1 ^ (u) = (u) = 1 The inset shows an enlargement of the plot. 



We express all quantities of interest in terms of the customary dimensionless (i.e., "re- 
duced") units. For example, length is given in units of a, energy in units of e, and temper- 
ature in units of e/k-Q. Other derived quantities are expressed in terms of suitable combi- 
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nations of these basic quantities. For example, stress is given in units of e/a 3 . Throughout 
this work we fix the temperature T = 1.0 which should be sufficiently subcritical taking the 
mean- field phase diagram of Hess and Su as a rough guidance [^(J. Our systems comprise 
N = 3000 fluid molecules which we equilibrate during an intial 10 4 Monte Carlo cycles. To 
save computer time we cut off fluid-fluid interactions beyond a separation of r c = 3.0 be- 
tween the centers-of-mass of a pair of fluid molecules. In addition, we employ a linked-cell in 
combination with a conventional (Verlet) neighbor list as described in the book by Allen and 
Tildesley to further speed up the simulations (see Chap. 5.3 of Ref. |22|). This list includes 
all particles as neighbors whose centers-of-mass are located within a distance of r n = 3.5 
from a reference molecule at the origin of the neighbor sphere. 

The simulations of this work are carried out for a substrate separation s z0 = 20. This value 
has been chosen because then our confined fluid contains a sufficiently large region centered 
around the middle of the slit-pore in which fluid molecules do not interact with either of the 
two substrates. This is illustrated by plots in Fig. 13. II indicating that (z, u) (z, u) < 
3 ■ 10~ 3 for \z/s z0 \ < 0.20. 



B. Properties 



To investigate the impact of the various hybrid anchoring scenarios introduced in Sec. IHBI 
we introduce in the following some key properties on which we shall base our subsequent 
discussion. The most important one is the so-called alignment tensor 23 1 



1 N 

1=1 



Ui 



(3.1) 



where "(g)" denotes the direct (i.e., dyadic) product and 1 is the unit tensor. The alignment 
tensor can be represented by a real, symmetric, traceless, 3x3 matrix that can be diago- 



24| . For a system with biaxial symmetry it can 



nalized numerically using Jacobi's method 
be shown that the diagonalized alignment tensor can be written as [see eqn. (lAllj) 

( -A m /2-C ^ 
-A m /2 + C 
A r 



diag Q = A 



V 



(3.2) 



/ 



This expression agrees with the one derived 



matrices 



25]. Following previous workers 



w Low who uses an expansion in terms of Wigner 
7, Q, Q we adopt S = (A m ) as a definition of 



S 



the Maier-Saupe nematic order parameter 28|, |29J] where the angular brackets denote an 



ensemble average in the specialized isostress ensemble. Once the largest eigenvalue of Q is 
known one may compute £ from eqn. (13.21) and with it the biaxial order parameter £ = (£). 

Other structural quantities that we shall be considering below take notice of the 
anisotropy and inhomogeneity of liquid crystals in confinement. The simplest one of these 
is the local density of the confined fluid defined as 

\ 

/ 6z \ A 

where 5 (z — Zj) is the Dirac 5-function and N(z) is the number of molecules within a small 
interval 5z centered on z. In addition, we consider the local nematic order parameter S(z) 
which we compute from an expression for Q (z) analogous to eqn. ( 13.11) replacing, however, Ui 
and N by their local counterparts Ui(z) and N(z), respectively. However, to gain even deeper 
insight into the nature of nematic phases it turns out to be useful to consider another measure 
of local nematic order provided by the (local value of the) second Legendre polynomial 

1 / 1 N{z) \ 
^W^2\ivM^ [3( ^' gz)2_1 V (3 ' 4) 

Hence, P z (z) = 1 if all molecules are homeotropically aligned (\v,i • e z \ = 1). If, on the 
other hand, molecules are aligned in a parallel fashion P z (z) = —0.5 (|itj • e z | = 0). Notice, 
that in writing eqn. ( 13. 4 j) we assume n to point along the z-axis of a space-fixed Cartesian 
coordinate system whereas the true direction of the local director n (z) computed as an 
eigenvector of Q (z) is unrestricted. 

In addition to these structural quantities we compute a specialized isostress heat capacity 

5, (H 2 )-(U) 2 
C - = 2 fcB+ Nk B T 2 (3 ' 5) 

that we already introduced in our previous study [7j. As we have demonstrated there, c T 
is particularly useful to identify characteristic stresses at which the confined liquid crystal 
may undergo subtle structural transformations. In eqn. (13. 5p 

H(R, U) = U(R, U) - t\\As zQ (3.6) 

such that 1(H) is a specialized enthalpy. 
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C. Heat capacity and global order parameters 



We begin our discussion with plots of c T as a function of applied compressional stress and 
various homogeneous and hybrid anchoring scenarios in Fig. 13.21 The isotress heat capacity 
exhibits pronounced peaks at sufficiently large compressional stress. As we have already 
demonstrated in our previous work, maxima of c T are fingerprints of an IN transition in 
the confined phase. Therefore, the stresses at which these maxima are located depend on 
the specific homogeneous anchoring scenario. For example, the plot in Fig. 13.2( a) shows 
that the IN transition occurs at lower compressional stress |rJ ~ 1.50 for the homogeneous 
homeotropic compared with Ty — 1.60 for the homogeneous planar anchoring scenario. 
Moreover, these IN transition occurs at a somewhat lower stress than Jry | ~ 1.70 observed in 
the bulk under identical thermodynamic conditions (see, for example, Fig. 4.3(a) of Ref. 3). 
This reflects the supportive character of the solid substrate in the formation of nematic 
phases. Note that in our previous study {t| we have been unable to observe an IN transition 
for the (homogeneous) homeotropic anchoring scenario because at the smaller pore width 
s z o = 10 considered there the impact of the walls is so strong that even at vanishingly 
small compression almost the entire confined fluid is already in a nematic phase. As we also 
explained earlier the height of peaks in the plots of c T reflect a loss of rotational entropy at 
the IN transition which is different for different anchoring scenarios (see Sec. V of Ref. |7J). 

Comparing results for the two homogeneous anchoring scenarios in Fig. 13.2( a) with the 
corresponding hybrid one two main differences can be recognized. First, the dominant peak 
in c T is shifted to a compressional stress rJ — 1.75 exceeding that characteristic of the 
two homogeneous anchoring scenarios. Obviously, the heterogeneity of the anchoring at the 
two substrates inhibits the formation of nematic phases. This seems quite striking at first 
in view of the fact that the direct impact of the solid substrates is restricted to a distance 
of about Az ~ 4 from either substrate surface as plots in Fig. 13.11 suggest. However, one 
has to keep in mind that nematic phases are characterized by long-range rotational order. 
Molecules of the confined soft-matter phase are therefore cooperative in transmitting the 
specific orientation imprinted by a solid substrate over distances substantially exceeding 
the range of the fluid-substrate interaction itself. The inhibition of the IN transition in 
the hybrid anchoring scenario may therefore be interpreted as frustration of the molecules 
in the confined phase which are trying to adjust their orientation with respect to the two 
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FIG. 3.2: Isostress heat capacity c T as a function of compressional stress ry for various homogeneous 
(□, O) and corresponding hybrid (•) anchoring scenarios, (a) Homeotropic (□), planar (O). (b) 
Nonspecific (□), planar (O). (c) Nonspecific (□), directional (O). Insets are enlargements of plots 
for hybrid anchoring scenarios. Solid lines are fits intended to guide the eye. 
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FIG. 3.3: "Snapshot" of a representative configuration in the nematic phase. Gray boxes indicate 
the position of the substrates in the x-z plane where the specific anchoring scenario is indicated 
within each box. The purple lines indicate the (virtual) boundaries of the simulation cell at which 
periodic boundary conditions are applied. To set the color code the global director is arbitrarily 
taken to be identical with the substrate normal. The sequence blue — > green — > red indicates a 
change in molecular orientation from parallel to perpendicular with n . 

conflicting anchoring conditions simultaneously. The structural frustration is illustrated in 
Fig. 13 .31 where we show a "snapshot" of a representative molecular configuration taken from 
the Monte Carlo simulations. 

A second interesting feature visible in the plot of c T in Fig. 13.2( a) for the hybrid anchoring 
is a smaller but significant secondary peak that occurs at lower compressional stress ry| — 
1.50 compared with the main peak at |ry| ~ 1.75. The existence of this secondary peak 
has been verified in Monte Carlo runs in which we enlarged the number of cycles up to 
3 x 10 6 cycles while carefully monitoring the evolution of c T as a function of the number of 
cycles similar to plots presented in Fig. 4.2 of Ref. Q- Comparing the plots for hybrid and 
homogeneous homeotropic alignment in Fig. 13.2( a) reveals that the secondary peak in the 
plots of c T for the hybrid anchoring coincides with that in the corresponding curve for the 
homogeneous homeotropic alignment. The origin of the secondary peak will become clear 
in Sec. IIIIDI where we discuss local structural features of the confined liquid crystals. 

For the hp anchoring scenario we also investigated the local density p (z) [see eqn. ( 13. 3p ] 
for a state point in the nematic phase. The plot in Fig. 13.41 shows that in the confined liquid 
crystal molecules arrange themselves in individual layers in the vicinity of either substrate 
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FIG. 3.4: Plot of local density p (z) as function of position between homeotropically (z w /s z = —0.5) 
and planar anchoring substrate (z w /s z = +0.5) at |tii| = 2.00. The thin solid line is intended to 
guide the eye whereas the horizontal thick solid line corresponds to ~p in the bulk-like region of the 
confined phase. The inset shows a plot of Ap where straight lines are fits of the far right side of 
eqn. (|3.T() to p(z m ) (see text) near the homeotropically (•) and planar (O) anchoring substrate. 



indicated by the oscillatory character of p (z) in the vicinity of the substrates. As one moves 
towards the midplane of the slit-pore the fluid becomes homogeneous (i.e., bulk-like) such 
that lim p (z) ~ p ~ 0.88 is independent of z. As long as the confined liquid crystal is fluidic 

z->0 

one expects 



Ap ee p (z m ) - p = aexp (z/£ ) cos (&Z 



6>->n-/2 
«1-V0 



aexp (z/£ ) 



(3.7) 



where z m is the position of maxima and p is the mean density of the bulk-like region in the 
plot of p (z). The correlation lengths £o and £i are determined by bulk prop erties whereas a 



and 9 depend on the specific nature of the fluid-substrate interactions |30| . Note, however, 
that eqn. ( 13. 7p holds strictly only for spherically symmetric molecules but is assumed here 
to be approximately valid because our system is composed of weakly anisometric molecules. 
A fit of eqn. ( 13. TP to p(z m ) gives £i ~ and 9 ~ n/2 such that the far right side of the 
equation provides an excellent representation of our data (see inset in Fig. I3.4p . Moreover, 
the fit yields £q ~ 1.5859 near the homeotropically anchoring substrate and £q ~ 1.2544 
near the planar anchoring one. The ratio £q/£o — 1-26 is nearly identical with the aspect 
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ratio of the ellipsoidal fluid molecules as one would have anticipated. 

As we already explained above the homogeneous nonspecific anchoring scenario favors 
a homeotropic alignment of the molecules in the nematic phase as well. However, as one 
can see from the corresponding plots in Fig. 13.2( b) the maximum in the plot of c T for this 
anchoring scenario is shifted to higher compressional stress ~ 1.60 compared with its 
counterpart for homogeneous homeotropic anchoring in Fig. 13.2( a). This indicates that the 
substrates for the nonspecific homogeneous anchoring scenario are somewhat less supportive 
in the formation of nematic phases. In fact, in the present case maxima of c r for the two 
homogeneous (nonspecific and planar) anchoring scenarios are located at about the same 
compressional stress as the plots in Fig. 13.2( b) indicate. 

The plot of c T for the np anchoring scenario exhibits qualitatively the same features 
already observed in Fig. 13.2( a) for the hp anchoring scenario. That is c T exhibits a main 
maximum at a larger compressional stress |t||| ~ 1.70 and a smaller secondary maximum 
at about rJ ~ 1.63. The shift of the secondary peak from ItJ ~ 1.50 [hp anchoring, 
see Fig. 13.2( a)] to \t\\\ ~ 1.63 [np anchoring, see Fig. 13.2( b)] is more pronounced than the 
shift of the associated main peak in the heat-capacity curves. Notice also that the location 
of the secondary maximum in the plot of c T for the two hybrid anchoring scenarios in 
Figs. 13.2( a) and 13.2( b) correlates nicely with the values \t\\ \ ~ 1.50 and my — 1.60 for the 
homogeneous homeotropic and nonspecific anchoring scenarios, respectively. This suggests 
that the secondary peak visible in the plots of c T for the two hybrid anchoring scenarios in 
Figs. 13.2( a) and 13.2( b) reflects an increase in homeotropic orientation of the molecules prior 
to the IN transition. This notion is supported by a more detailed analysis of local structural 
features presented below in Sec. IIIID1 

However, before turning to that discussion it is instructive to consider yet another hybrid 
anchoring scenario which combines nonspecific anchoring at one substrate with directional 
anchoring at the other one. As plots in Fig. 13.2( c) show the maximum in c T for the homo- 
geneous directional anchoring occurs at a compressional stress of about |t||| ~ 1.45 which is 
the smallest value of the three homogeneous anchoring scenarios considered in this study. 
This indicates that directional (planar) anchoring is the most supportive scenario in the for- 
mation of nematic phases in confined liquid crystals. Combining directional with nonspecific 
anchoring causes the small but pronounced secondary peak in the plots of c T in Figs. 13.2( a) 
and 13.2( b) to give way to a broad shoulder at about \t\\\ ~ 1.55 which is located approxi- 
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FIG. 3.5: The Maier-Saupe nematic order parameter S as a function of applied compressional 
stress. Homogeneous anchoring scenarios: (x) planar, (+) nonspecific, (*) directional; hybrid 
anchoring scenarios: (□) np, (O) nd. Also shown are plots of the biaxial order parameter £ for np 
(■) and nd anchoring scenarios (•). Solid and dashed lines represent fits to guide the eye. The 
vertical dashed line indicates typical values of £ obtained for homogeneous anchoring scenarios. 

mately halfway in between the peaks of c T at rj| — 1.45 for the homogeneous directional 
and |t||| ~ 1.60 for the homogeneous nonspecific anchoring scenarios. One may therefore 
speculate that the broad shoulder in the plot of c T in Fig. 13.2( c) reflects some other struc- 
tural transformation than the secondary peaks in the plots of c T in Figs. 13.2( a) and 13.2( b) 
(see Sec. HUE) . 

The IN transition may also be located through plots of the (global) Maier-Saupe nematic 
order parameter S plotted as a function of applied compressional stress in Fig. 13.51 In all 
cases considered S(t\\) is sigmoidal in shape. Taking the inflection of the sigmoidal curves as 
an operational definition of the location of the IN transition j3] very good agreement with 
the position of maxima in the plots of c T versus r\\ plotted in Fig. 13.21 is obtained. For the 
two hybrid anchoring scenarios the curves S(t\\) are indistinguishable which is consistent 
with the observation that the main peak of c T in the parallel plots in Figs. 13.2( b) and 13.2( c) 
appears to be unaffected by the specific combination of anchoring conditions at the two 
substrates. In general, it turns out that in the nematic phase S is much smaller for hybrid 
than for homogeneous anchoring scenarios. However, the specific hybrid anchoring scenario 
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manifests itself in plots of the biaxial order parameter £ which turns out to be larger for 
the nd than for the np anchoring scenario. Also shown in Fig. 13.51 are typical values of £ for 
homogeneous anchoring where one expects £ = 0. However, in the actual simulations one 
typically obtains £ < 0.02 on account of a small system-size effect discussed in detail in the 



Appendix of the paper by Eppenga and Frenkel 261 ] . 



D. Local orientation 

To gain deeper insight into structural transformations occuring in the confined liquid 
crystal with increasing compressional stress we consider the local Maier-Saupe nematic or- 
der parameter S (z) and Pi {£) [see eqn. (13.41) ]. According to their definition P 2 (z) = S (z) if 
molecules are aligned homeotropically with the solid substrate. If this alignment is perfect, 
P2 (z) = 1- An inspection of Fig. 13.61 shows that this is the case in the immediate vicinity 
of the lower substrate if g (u) is given by eqn. ( I2.8bl) . If, on the other hand, the preferential 
orientation of molecules is in a direction parallel with the solid substrates, P2 {z) = —0.5 but 
S (z) > which is the case in the immediate vicinity of the upper substrate characterized 
by the anchoring function specified in eqn. fl2.8c[) . For the plots presented in Fig. 13.61 we se- 
lected state points around the two maxima visible in the plots of c T in Fig. 13.2( a) for hybrid 
anchoring. Both parts of Fig. 13.61 show that the system undergoes a structural transition 
during which the local order increases. These transformations are restricted mostly to "in- 
ner" parts of the confined fluid (i.e., regions located at \z\ /s z q < 0.25) for which the direct 
interaction with either solid substrate is negligible according the plot in Fig. 13.11 Notice 
also that in the "nematic" phase the local order parameter S (z) is still smaller than the 
threshold value of 0.4 suggested by the Maier-Saupe (mean-field) theory for the IN phase 
transition for z/s z > 0.1 [see Fig. 13.6( b)]. This is the reason why the global nematic order 
parameter S for the two hybrid anchoring scenarios plotted in Fig. 13.51 barely reaches the 
value S — 0.4 even at the highest compressional stresses considered. 

To visualize subtle details of these structural transformations more clearly it turns out 
to be helpful to consider AS (z) = S(z; T|j 2 ^) — S(z; r,| ) and, with an analogous definition, 
AP 2 (z) where |r,| | < l r jj 2 ^|- Plots of both quantities in Fig. 13.7( a) show that around the sec- 
ondary maximum of the heat capacity curve presented in Fig. 13.2( a) for the hybrid anchoring 
scenario AS (z) = P2 (z) with a maximum located closer to the homeotropically anchoring 
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FIG. 3.6: Plots of S (z) (lines) and P2 (z) (symbols) as functions of position between lower and 

upper substrates located at z w /s z o = =F0-5 (s z o = 20) for the hp anchoring scenario, (a) ( ), (A): 

Tii = -1.45; (---), (O): T » = -1.55. (b) ( ), (A): tn = -1.70; (---), (O): r,, = -1.80 [see also 



Fig. 13.2( a)], The full horizontal line demarcates S (z) = 0.4 taken as the nematic threshold 



substrate. For thermodynamic states around the main maximum of c r in Fig. 13.2( a). on 
the other hand, plots in Fig. 13.7( b) indicate that the maximum in the curves AS" (z) and 
AP 2 (z) is shifted towards the parallel anchoring substrate whereas AS (z) = AP 2 (z) still 
holds. 

To interpret the data presented in Fig. 13 .71 it is useful to consider the schematic represen- 
tation presented in Fig. 13.81 In the nematic phase the director n describing the preferential 
orientation of the molecules will generally form a certain angle with the substrate normal 
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FIG. 3.7: Plots of AS* (z) (•) and AP2 (z) (O) for the hybrid homeotropic-planar anchoring scenario 
as functions of position between lower and upper substrates, (a) tP = -1.45, rj 2) = -1.55, (b) 
Ty 1 ^ = —1.70, T|j 2 ^ = —1.80. Notice, that error bars may increase substantially as \z\ /s z o — > 0.5 
because of the excluded volume in the immediate vicinity of the solid substrates (see Fig. 13. 6p , 

(i.e., the z-direction) and may lie anywhere on the blue cone drawn in Fig. 13.81 The projec- 
tion n z of n onto the z-axis is related to the value of P2 according to the definition given in 
eqn. f)3.4p . Because \n\ = 1 by definition, the radius R of the circular base area of the blue 
cone becomes smaller and vanishes eventually if 9 — > 0. If, on the other hand, n becomes 
more aligned with the x-y plane, R increases and n z vanishes if ip — > 0. Therefore, R may 
be viewed as a geometrical representation of the degree of homeotropic alignment. Likewise, 
the radius r of the circular base area of the yellow cone is a geometrical representation of 
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FIG. 3.8: The director n, which ends at the center of the base area of the yellow cone, in the 
nematic phase forms an angle cos 9 = n • e z = ra z with the z-axis, an angle <p = ^ — 6 with the 
x—y plane, and may lie anywhere on the surface of the blue cone. Radius R of the base area of the 
blue cone is a measure of the degree of homeotropic alignment proportional to P 2 [see eqn. (|3.4p ]: 
radius r of the base area of the yellow cone is a measure of S. 

the Maier-Saupe nematic order parameter S: The larger r the lower is S. A nonvanishing 
value of r reflects a certain width of the distribution of molecular orientations around n on 
account of thermal fluctuations. 

Applying these geometrical concepts to the plots in Fig. 13.7( a) it is clear that R and r 
of the two cones plotted in Fig. 13.81 shrink simultaneously as the transverse compressional 
stress increases. This implies that the nematic order parameter S (z) increases and that the 
director becomes more aligned with the z-axis. In other words, the secondary maximum in 
the plot of c r in Fig. 13.2( a) at |t||| ~ 1.50 reflects an increasing homeotropic (pre) alignment 
of molecules in an isotropic region of the confined fluid [cf., Figs. 13.6( a). 13.7( a)], This 
homeotropic prealignment takes place in the part of the fluid closer to the homeotropically 
anchoring substrate. Similarly, the maximum at c r at Ty ~ 1.75 may be interpreted as an 
increase in nematic order which is mostly caused by an increasing alignment of molecular 
orientation with the z-axis. However, here we note that in the vicinity of the planar anchoring 
substrate AS (z) < AP 2 (z). In other words, the increasing homeotropic alignment near the 
planar anchoring wall is associated with a slight loss of nematic order in this region. 
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FIG. 3.9: As Fig. 13.71 but for the np anchoring scenario, (a) = —1.58, T^p = —1.65; (b) 
tP = -1.67, r{ 2) = -1.80. 

Corresponding plots in Fig. 13.91 for the np anchoring show that around the secondary 
peak of c T at rn ~ —1.63 [see Fig. 13.2( b)] the structural transformation is very similar 
to the one illustrated by the plots in Fig. 13.7( a): molecules prealign preferentially with 
the z-axis in that region of the confined fluid controlled by the nonspecifically anchoring 
substrate. In the geometrical interpretation of Fig. 13.81 this may again be desribed in terms 
of a simultaneous shrinkage of both R and r. However, the plots in Fig. 13.9( b) illustrate a 
somewhat different structural transformation. Here, AS (z) > AP 2 {z) indicating that the 
alignment of molecules with the z-axis is smaller than the increase in nematic order. In the 
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geometrical terms of Fig. 13.81 this may be desribed as a decrease of R and a simultaneous 
larger decrease of r such that 9 decreases slightly in an increasingly nematic phase. 
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FIG. 3.10: As Fig. 13.71 but for the nd anchoring scenario, (a) = —1.50, 
rj 1 ) = -1.63, r,( 2) = -1.75. 



■1.60; (b) 



A totally different structural transformation is observed for the nd anchoring scenario for 
which we present plots of AS* (z) and AP 2 ( z ) i n Fig. 13.101 For two thermodynamic states 
in the vicinity of the weak shoulder in the plot of c T at tJ ~ 1.55 [see Fig. 13.2( c)] the plots 
in Fig. 13.10( a) indicate that in the immediate vicinity of the hybrid anchoring substrate 
AS (z) = AP2 (z) referring to an enahnced homeotropic alignment of the molecules as the 
compressional stress increases. On the contrary, in the upper halfspace of the slit-pore (i.e., 
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for z/s z q > 0), AP2 (z) < whereas AS (z) > 0. In the geometric terms introduced in 
Fig. 13.81 this can be ascribed to an increase of R (i.e., a decrease of n z ) and a decrease of r. 
In other words, the increase in nematic order is a consequence of a preferential alignment of 
the molecules with the x-axis. At higher compressional stresses in the vicinity of the main 
peak of c r at |t|| | ~ 1.70 [see Fig. 13.2( c)] this trend is somewhat reversed as the corresponding 
plots in Fig. 13.10( b) show. The plot of AS" (z) is symmetric about the midplane of the slit- 
pore and because AS* (z) > the nematic order in the system increases mainly at the center 
of the confined fluid. There is also some enhanced homeotropic alignment which turns out 
to be somewhat more pronounced in the vicinity of the homeotropically anchoring wall in 
accord with one's physical intuition. 

IV. SUMMARY AND CONCLUSIONS 

In this work we consider a simple model of a liquid crystal confined to a slit-pore with 
hybrid anchoring at the solid substrates, that is each substrate discriminates molecules on 
account of their orientation with respect to the substrate plane where we take this energetic 
discrimination to be different at the two substrates. Employing Monte Carlo simulations in 
a specialized isostress ensemble we can locate specific compressional stresses, at which the 
fluid undergoes a significant structural reorganization, through characteristic maxima of the 
isostress heat capacity. 

As in our previous study |3] the location of these maxima coincides with inflection points 
in plots of the Maier-Saupe nematic S order parameter as a function of compressional 
stress. Considering various hybrid anchoring scenarios it turns out that the location of 
the IN transition is largely independent of the combination of anchoring conditions at the 
two substrates. This is because in the nematic phase the director is trying to point in a 
direction compromising between the competing orientations enforced by the two substrates 
simultaneously. Consequently, plots of S versus r\\ are nearly indistinguishable as far as the 
np and nd anchoring scenarios are concerned. Our parallel analysis of local order parameter 
profiles shows that the IN transformation occurs near the center of the confined liquid crystal 
that is approximately halfway in between the solid substrates. 

However, the specific hybrid anchoring scenario does have an impact on a secondary 
structural reorganization that occurs at lower compressional stresses preceding that of the 
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IN transition. This secondary transformation may cause prealignment with the z-axis (i.e., 
in a direction perpendicular to the substrate plane) or it may force molecules to prealign to 
a certain extent with the x-y plane. Which one of these two cases is realized depends on 
the stronger anchoring scenario in the hybrid pair. Thus, our results show that homeotropic 
is stronger than planar anchoring during this prealignment process. The same is true for 
the nonspecific anchoring which favors a homeotropic alignment of the fluid molecules on 
account of the symmetry breaking presence of the solid substrates. In the combination with 
nonspecific anchoring the directionally anchoring substrate "wins" and is thus considered 
stronger in the above sense. 

In slit-pores with hybrid anchoring at the substrates one generally anticipates a certain 
degree of biaxiality. Our simulations do, however, show that the degree of biaxiality is 
actually quite small. Only the nd anchoring scenario exhibits a significant degree of biaxiality 
of about 25% of S in the nematic phase. The apparent lack of a more substantial biaxiality 
can be understood by realizing that the director field n (z) varies as a function of position 
w.r.t. the location of the substrate plane. This is inferred from a characteristic "snapshot" 
of a configuration of molecules in the nematic phase of the nd anchoring scenario. In other 
words, the preferential orientation of molecules in the nematic phase varies with position 
such that more than just two preferred orientations in the confined liquid crystal are present. 
The small nonvanishing biaxiality for the nd anchoring scenario is caused by those few layers 
of molecules in the vicinity of each substrate that exhibit either a homeotropic or a planar 
orientation of the molecules. 

Appendix A: Perturbational treatment of biaxial symmetry 

In this appendix we derive eqn. (13.21) by perturbational arguments. We assume that we 
can split the alignment tensor Q for a system with biaxial symmetry into a contribution 
q(o) representing a system with uniaxial symmetry (i.e., the "symmetry ground state") and 
a perturbation which removes the rotational symmetry in the plane orthogonal to the 
symmetry axis in the uniaxial case ("symmetry excited state"). Hence, we may write 

Q = Q(°) + Q« (Al) 
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where Q, Q^, and are assumed to be represented by symmetric, traceless, 3x3 
matrices. Let satisfy the eigenvalue equation 



(o)-(o) 



Q w x 



A (o)^(o) 



1,...,3 



(A2) 



where {ccf^} and {A 4 W } are the sets of eigenvectors and -values of Q^, respectively where 



,(o) 



we introduce matrix elements of Q w via Q\j' = xf'Q^x^'. Because Q w is real and 
symmetric, Qfj 



(o) 



;(°)n(0)^(°) 



Qji . It then follows directly from eqn. (|A2[) that 



A (o) 



X; 



(0) 



(A3) 



which expresses the well-known fact that eigenvectors to different eigenvalues are automat- 
ically orthogonal to one another. Let X^ = A m be the largest eigenvalue of Q^. Uniaxial 



symmetry then requires A^ 

;(0) 



(0) 



A 



(o) 



A' < A m implying x ^ ■ xf' 







1, 2), that is 



eigenvector x 3 is orthogonal to the other two. However, because A^ 



(0) 



A 



(0) -(0) 



. X cLlld Xr) 



(0) 



can lie anywhere in the plane orthogonal to xf\ In particular, eqn. ( 1A3I) remains valid 
without requiring these two eigenvectors to be orthogonal to one another. Hence, only x^ 
defines a distinct direction and represents the axis of uniaxial symmetry. 

To proceed it is important to realize that as long as the original eigenvectors are linearly 
independent one may apply a standard orthogonalization procedure (such as, for example, 

Sij, Wi,j. 



Gram-Schmidt, see Ref. 



32l ) to replace {$,-°' ) } by {y^} where now yf ■ y 



(o) 



Moreover, members of the new set {y^} are still eigenvectors of to the same eigen- 
values A' and A m . Because {y^} form a complete orthonormal set of eigenvectors one may 
construct a matrix Y where column i is formed by yf 1 ^- With this matrix Q(°) may be 
diagonalized according to 



( X' ^ 



diagQ(°) = Ao = Y - 1 Q(°)Y 



V 



A' 
A m 



/ 



/ 



V 



-A m /2 
-A m /2 
A„ 



/ 



(A4) 



where the far right side follows because the trace is an invariant of a matrix such that 
TrQ(°) = TrA = 0. 

Assume now that Q satisfies an eigenvalue equation similar to eqn. ( 1A2I) that is 



Qy = Xy 



(A5) 
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where we may express the eigenvector y in the basis {y\ } according to 

^ ^(0) . ^(o) . -(0) /Ar x 
y = a x y\ + a 2 y\ + a 3 y y 3 (A6) 

We determine the expansion coefficients {aj} by inserting eqn. (IA6I) into eqn. (IA5[) and 
multiply the resulting expression from the left successively by yf \ i — 1, . . . , 3. Employing 
the orthonormality of the basis functions we obtain a set of three coupled linear equations 
which has a nontrivial solution if and only if the secular determinant 

-A m /2 - A ±C 
±C -A m /2 - A 
A m - A 

The specific form of the secular determinant follows from eqn. (lAip and (1A2I) assuming that 



(A7) 



My 



o 



(0) 



(A8a) 
(A8b) 



where £ is the biaxiality parameter. We have chosen the perturbation such that it preserves 
the orthonormality of the basis function but removes the degeneracy of the eigenstates in 
the plane orthogonal to xf^ = y^ by introducing a second length scale represented by (. 
This can be seen more directly by computing the roots of the third-order polynomial in A 
represented by eqn. (1A7I) which turn out to be given by 



A 
A± 



A 



in 
A r 



2 s 



where the associated eigenvectors of Q are given by 



2/o = 2/ 3 



(0) 



y ± = -^(y<°>±Sf) 



(A9a) 
(A9b) 

(AlOa) 
(AlOb) 



Using the eigenvectors yo and y± as a new basis we may set up a matrix Y analogous to 
Y in eqn. (1A4|) and diagonalize Q according to 



/ 



diagQ = Y^QY 



-A m /2 + C 
-A m /2-C 
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(AH) 
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